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Abstract 

This  paper  documents  the  investigation  of  fuel  injector  fault 
detection  methods  for  a  seven  liter  diesel  engine.  This  effort 
was  conducted  for  the  Tank  Automotive  Research 
Development  Engineering  Center  (TARDEC)  Condition 
Based  Maintenance  (CBM)  team.  The  task  was  to  develop 
algorithms  capable  of  real-time  detection  of  injector  misfire 
events.  The  purpose  of  this  task  was  to  enable  TARDEC ’s 
Engine  Control  Management  (ECM)  research  and 
development  efforts  to  evaluate  the  technical  feasibility  of 
integrating  automated  on-board  condition  monitoring 
algorithms  with  future  ECM  monitoring  and  control 
operations. 
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During  this  investigation,  it  was  shown  that  multiple 
techniques  can  correctly  detect  and  identify  injector  cylinder 
misfiring.  Each  individual  technique  has  its  own 

advantages,  and  this  investigation  focused  on  low 
computational  power  signal  processing  methods  that  would 
be  suitable  for  embedding  in  an  engine  controller  or 
processor.  The  list  below  includes  six  injector  fault  analysis 
approaches  that  were  evaluated  for  this  effort: 

1.  Injector  signal  based  analysis 

2.  Vibration  based  analysis 

3.  Time  domain  speed  encoder  analysis 

4.  Order  domain  speed  encoder  analysis 

5.  FFT  classifier  selection  techniques 

6.  Time  domain  classification  techniques 

This  paper  highlights  the  first  four  diagnostics  techniques. 
On-platform  tests  are  suggested  for  technique  validation  and 
future  development  of  these  initial  findings. 

1.  Introduction 

The  Tank  Automotive  Research  Development  Engineering 
Center  (TARDEC)  Condition  Based  Maintenance  (CBM) 
team  investigated  the  development  of  a  Vehicle  Health 
Management  Systems  (VHMS)  for  ground  combat 
platfoms.  The  mission  of  the  program  was  to  direct  the 
development  of  technologies  to  provide  the  US  Army  with 
an  improved  diagnostic,  predictive  and  sustainment 
capability  for  ground  vehicles. 
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In  2010,  TARDEC  conducted  diesel  engine  dynamometer 
testing  on  a  seven  liter  diesel  engine.  As  part  of  the 
dynamometer  test  schedule,  an  engine  testing  tool  was  used 
to  selectively  deactivate  individual  fuel  injectors  over  a 
range  of  speed  and  load  conditions.  This  set  of  tests  formed 
the  basis  for  the  data  set  used  in  this  investigation. 

This  effort  investigated  multiple  approaches  and 
documented  the  advantages  and  limitations  with  varying 
degrees  of  robustness,  cost  and  complexity  in  terms  of 
instrumentation,  hardware  and  software  requirements.  From 
failure  mode,  effects  and  criticality  analysis  (FMECA) 
studies  of  diesel  engines  shown  by  Banks,  Hines,  Febold, 
Campbell,  Begg  and  Byington  (2001)  can  be  used  to  detect 
failures  that  will  produced  the  greatest  value  for  achieving 
high  operational  availability  and  reducing  maintenance 
costs.  For  this  effort,  we  focused  on  detecting  injector  (or 
combustion)  related  faults.  We  emphasize  the  distinction 
that  unlike  the  work  developed  in  the  past  by  Murphy, 
Febold,  Reichard,  Galie,  and  Byington  (2003)  the 
approaches  described  here  focused  on  commonly  installed 
sensors.  Moreover,  the  goal  of  this  work  was  to 
developsimple  data  collection/processing  schemes  for 
embedded  solutions. 

The  six  injector  fault  analysis  techniques  studied  as  part  of 
this  effort  can  be  grouped  into  three  domains:  time, 
frequency,  and  classification.  Each  technique  is 
summarized  below. 

Time  Domain  Techniques: 

■  Injector  Signal  Based  Analysis 

o  Begin  with  calculating  the  timed  based  injector 
pulse  profile. 

o  Once  the  pulse  profile  is  characterized,  determine 
injector  faults  by  categorizing  pulse  profile  energy. 

■  Vibration  Based  Analysis 

o  Investigate  an  alternative  accelerometer-based 

approach  utilizing  only  analog  acquired  signals, 
o  Similar  to  the  injector  signal  analysis,  this 

approach  also  determines  injector  faults  by  the 
presence  of  energy  in  the  vibration  signal  during 
the  specific  cylinder  firing  event  window. 

•  Time  Domain  Speed  Encoder  Analysis 

o  Investigate  fault  detection  approaches  using  the 
crankshaft  speed  encoder  signal  acquired  from  the 
digital  timer  board  DAQ  system, 
o  Determine  fault  based  on  the  instantaneous  speed 
profile  of  the  crank  shaft. 


Order  Domain  Techniques: 

•  Order  Domain  Speed  Encoder  Analsis 

o  Investigate  the  order  domain  analysis  of  the 
crankshaft  speed  encoder  signal  acquired  from  the 
timer  board  DAQ  system. 

o  Determine  fault  based  on  the  magnitude  and  phase 
at  particular  crankshaft  orders. 

Classification  Domain  Techniques: 

•  General  Fault  Classification  and  Reasoning 
Algorithms 

o  Investigate  fault  detection  methods  using 

classification/reasoning  algorithms  on  FFT 
frequency  bin  data. 

•  Time  Based  Classification  Approach 

o  Investigate  fault  detection  methods  using 

classification/reasoning  algorithms  on  time  domain 
data. 

For  this  paper,  only  the  time  and  order  domain  techniques 
will  be  discussed.  The  following  sections  will  discuss  the 
process  for  calculating  injector  signal  profile  and  the  four 
diagnostic  techniques  in  further  detail. 

2.  Engine  Injector  Analysis  Background 

The  purpose  of  injector  timing  analysis  is  to  quantify  the 
pulse  duration  and  delay  of  the  injection  of  fuel  into  the 
engine  with  respect  to  the  piston’s  orientation.  This  is 
usally  referenced  to  the  piston’s  Top  Dead  Center  (TDC) 
location  of  0  degrees.  Top  dead  center  is  the  position  of  a 
piston  in  which  it  is  farthest  from  the  crankshaft. 

In  graphic  terms  the  goal  of  this  task  is  shown  in  Figure  1. 
The  subsequent  sub-sections  describe  the  process  by  which 
this  end  result  is  obtained. 


Figure  1.  Graphic  Illustration  of  Fuel  Injector  Timing 
Analysis  Objective 


The  algorithm  developed  to  perform  this  task  calculates  this 
pulse  width  for  every  injector  actuation  event  that  occurs 
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during  the  entire  ten  second  recording.  The  algorithm  is 
designed  to  differentiate  and  properly  account  for  injector 
pulses  consisting  of  single,  double  or  triple  pulse  profiles. 

The  data  in  Figure  2  illustrates  single,  triple  and  double 
pulse  profiles  respectively.  Single  pulses  were  typically 
observed  at  relatively  low  speed  around  700  revolution  per 
minute  (RPM)  and  low  load  (30  lbs-ft)  conditions  during  the 
initial  warm-up  of  the  engine.  Triple  pulses  were  observed 
as  speed-load  conditions  were  increased.  Typically  the 
right-most  pulse  within  the  three  pulse  profile  would 
gradually  widen  as  speed-load  continued  to  increase  and 
eventually  ‘absorb’  the  center  pulse,  thereby  resulting  in  the 
two  pulse  profile  at  the  highest  speed  -  load  (2200  RPM  - 
560  lbs-ft)  condition. 


While  this  processing  technique  is  the  basis  for  the  first 
time-domain  analysis  approach,  it  is  also  useful  for  insuring 
that  the  engine  controller  is  performing  properly  and  that  the 
engine  is  load  balancing  across  the  engine  properly.  The 
following  section  discusses  the  first  time  domain  technique, 
which  parses  the  data  similarly  to  the  approach  just 
mentioned  and  looks  at  the  energy  within  each  combustion 
event  window. 

3.  Injector  Signal  Based  Analysis 

This  time-domain  technique  for  injector  fault  detection 
looks  at  the  presence  of  energy  in  the  injector  signal  during 
the  injection  firing  event  window.  Figure  3  illustrates  the 
method  by  which  signal  energy  is  determined  to  exist  in  the 
detection  window. 
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Figure  3.  Detecting  Injector  Signal  Energy  in  the  Window 
of  Interest 

Figure  4  is  a  graphic  representation  illustrating  how  the 
algorithm  identifies  a  cylinder  three  injector  misfire.  This 
plot  is  predicated  on  the  fact  the  algorithm  has  identified 
TDC  of  cylinder  one.  Once  the  cylinder  one  TDC  is  located, 
the  algorithm  measures  the  injector  signal  energy  in  each 
detection  window.  If  the  signal  energy  does  not  exceed  the 
threshold  condition,  the  algorithm  identifies  the  respective 
injector  as  a  misfire  fault. 


Figure  2.  (a):  Single  Injector  Pulse;  (b)  Triple  Pulse;  (c) 
Double  Pulse 

The  algorithm  stores  the  computed  pulse  duration  (AT)  and 
the  relative  delay/offset  from  Top  Dead  Center  (TDC)  for 
each  respective  pulse  profile  event  occurring  during  the 
record  set.  In  broad  terms  the  analysis  is  divided  into  three 
sub-tasks: 

1 .  Compute  detection  window 

2.  Determine  true  top  dead  center 

3.  Process  injector  pulses 


The  injector  signal  based  detection  approach  resulted  in 
100%  detection  rate  for  all  of  the  specific  load,  speed  and 
fault  conditions  tested.  No  false  positives  or  negatives  were 
generated  using  this  approach.  The  sensing  requirements 
for  this  approach  are  an  injector  signal  and  a  once  per 
revolution  signal  (1PPR/TDC  signal).  The  signal  processing 
approach  works  under  varying  conditions  regardless  of  the 
number  of  injector  pulses  produced  per  cylinder, 
(commonly  one,  two  or  three  pulses).  Furthermore,  due  to 
the  low  computation  processing  power  needed  for  this 
approach,  this  simplistic  technique  can  be  easily  embedded 
into  an  engine  controller  or  small  microprocessor. 
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The  next  method  applies  the  same  signal  processing 
methodology  to  accelerometer  (vibration)  data  collected 
from  a  location  mounted  on  top  of  the  engine  heads. 


Figure  4.  Detecting  a  Representative  Cylinder  Misfire  Event 


4.  Vibration  Based  Analysis 

The  motivation  for  investigating  a  vibration-based  fault 
detection  approach  is  to  provide  an  alternative  method  for 
detecting  combustion  related  faults.  Therefore,  such  an 
approach  could  be  used  for  engines  without  injector  firing 
signals  or  engines  having  mechanical  injectors.  Besides 
developing  a  technique  using  vibration  sensors,  the  task  also 
investigates  the  minimum  number  of  accelerometers  needed 
to  correctly  detect  a  fuel  injector  misfire  in  any  given 
cylinder.  Like  the  previous  detection  approach,  this 
investigation  was  conducted  on  an  engine  dynamometer  test 
stand  while  seeding  one  injector  fault  electronically. 

The  data  in  Figure  5  illustrates  a  representative  plot  over 
two  revolutions  of  the  crankshaft  at  an  operating  condition 
of  30  lbs-ft  @  700  RPM.  The  alternating  yellow  and  white 
regions  correspond  to  crankshaft  angle  of  revolution  during 
which  fuel  injectors  are  injecting  fuel  into  respective 
cylinders,  designated  in  firing  order  sequence:  #1-5 -3 -6-2-4. 


Raw  accelerometer  and  injector  signal  overlay 
Baseline  •  No  Fault;  Sequence  Nr.:  492;  Speed:  700  RPM-  Load:  30  LBS-FT 
Firing  order:  1-5-3-6-24 


Crank  angle  (Degrees  with  respect  to  start  of  Cylinder  1  window) 


Figure  5.  Accelerometer  and  Injector  Signals:  Baseline  -  No 
Fault  (30  lbs-ft@700  RPM) 

The  data  in  Figure  6  shows  a  representative  plot  with  an 
injector  misfire  event  in  cylinder  one  at  the  same  operating 
conditions  (30  lbs-ft  @  700  RPM).  As  a  note,  this  snapshot 
of  data  was  collected  during  the  start  of  engine  warm-up, 
hence  the  single  injector  profile. 


0  100  200  300  400  600  600  700 

Crank  angle  (Degrees  with  respect  to  start  oF  Cylinder  1  window) 


The  vertical  axis  overlays  the  injector  profiles  for  each  of 
the  six  cylinders  and  the  signal  acquired  from  the 
accelerometer  mounted  on  the  number  three  cylinder  head  . 
Note  the  clear  relationship  between  the  two  injector  pulses 
per  revolution  and  the  corresponding  two  pulses  detected  by 
the  accelerometer. 


Figure  6.  Accelerometer  and  Injector  Signals:  Fault 
Cylinder  One  (30  lbs-ft@700  RPM) 

Here  again,  accelerometer  three  is  the  channel  used  to  detect 
the  vibration.  This  plot  demonstrates  that  this  centrally 
positioned  accelerometer  is  capable  of  detecting  injector 
misfire  events  at  this  operating  condition  on  the 
dynamometer  test  stand. 
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The  data  in  Figure  7  illustrates  the  signal  obtained  from 
accelerometer  three  at  a  relatively  higher  speed  and  load 
operating  condition  (560  lbs-ft  @  2200  RPM).  The 
correlation  between  the  injector  pulses  and  vibration  signal 
is  not  as  apparent  as  it  was  for  the  30  lbs-ft  @  700  RPM 
case  but  with  simple  low-pass  filtering  the  firing  event 
faults  can  still  be  detected. 


Raw  accelerometer  and  injector  signal  overlay 
Baseline  -  No  Fault;  Sequence  Nr.:  495;  Speed:  2200  RPM-  Load:  560  LBS-FT 
Firing  order:  1-5-3-6-2-4 
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Figure  7.  Accelerometer  and  Injector  Signals:  Baseline-No 
Fault  (560  lbs-ft@2200  RPM) 

Nonetheless,  as  illustrated  in  Figure  8,  accelerometer  three 
with  low-pass  filtering  is  capable  of  detecting  a  cylinder  one 
injector  misfire  event.  Accelerometer  three  demonstrated 
similar  results  in  its  ability  to  detect  misfire  faults  in  the 
other  five  cylinders. 


Crank  angle  (Degrees  with  respect  to  start  of  Cylinder  1  window) 


The  data  in  Figure  9  demonstrates  that  the  signal  obtained 
from  accelerometer  three  is  capable  of  detecting  injector 
misfires  across  each  of  the  six  cylinders,  respectively.  As  is 
shown  in  Figures  5  -  8,  the  horizontal  axis  in  Figure  9 
indicates  crank  angle  revolution  over  two  complete 
revolutions.  The  vertical  axis  indicates  the  signal  obtained 
by  accelerometer  number  three  for  six  respective  faulted 
cylinders  in  firing  order,  as  read  from  top  to  bottom,  (#1-5- 
3 -6-2-4).  Because  the  vertical  axis  is  arranged  in  order  of 
firing  sequence  with  respect  to  crank  angle,  the  highlighted 
green  diagonal  box  creates  an  ‘identity  matrix’.  Stated 
differently,  using  injector  misfire  three  as  an  example,  by 
visually  tracing  an  imaginary  vertical  line  beginning  on  the 
horizontal  axis  at  approximately  300°  up  to  the  point  where 
the  line  intersects  the  diagonal  green  box,  the  reader  will 
observe  the  amplitude  of  the  accelerometer  signal  is 
noticeably  smaller  relative  to  the  preceding  (left  most)  and 
subsequent  (right  most)  portions  of  the  signal.  Tracing  an 
imaginary  horizontal  line  from  this  green  box  to  the  left 
until  intersecting  the  vertical  axis,  the  reader  will  observe 
the  vertical  axis  number  label  corresponding  to  the  cylinder 
that  experienced  the  injector  misfire. 
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Raw  accelerometer  #3  signal  overlay  for  6  separate  faulted  cylinder  events 
Fault  Sequence#:  457  462  458  460  461;  Speed:  1800  RPM-  Load:400  LBS-FT 
Firing  order:  1-5-3-6-2-4 


Crank  angle  (Degrees  with  respect  to  start  of  Cylinder  1  window) 


Figure  9.  Accelerometer  Three  Detection  of  Injector 
Misfires  Occurring  in  Each  of  the  Six  Cylinders:  (560  lbs- 
ft@2200  RPM) 

The  vibration  based  detection  approach  resulted  in  a  100% 
detection  rate  for  all  the  specific  load,  speed  and  fault 
conditions  tested.  The  signal  requirements  for  this  approach 
are  vibration  data  from  one  accelerometer  located  in  an 
optimum  location  and  a  once  per  revolution  signal  (1PPR  or 
TDC  signal).  This  signal  processing  approach  also  works 
under  varying  conditions  regardless  of  the  number  of 
injector  pulses  per  cylinder. 


Figure  8.  Accelerometer  and  Injector  Signal  Overlay:  Fault 
Cylinder  One  (560  lbs-ft@  2200  RPM) 


As  observed  in  the  plots,  the  approach  requires  only  one 
accelerometer  to  detect  all  cylinder  faults.  Having  stated 
this,  this  approach  would  have  to  be  validated  on  vehicles 
operating  during  extreme  conditions  to  determine  effects  of 
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externally  induced  vibrations.  If  validated,  this  simple 
approach  can  also  be  easily  embedded  into  an  engine 
controller  or  small  embedded  processor  based  system.  It  is 
worth  noting  that  different  engines  would  each  require  their 
own  respective  signal  amplitude  threshold  levels  in  order  to 
trigger  a  fault  detection  alert.  This  paper  demonstrates  the 
feasibility  of  the  approach.  Amplitude  thresholds  were  not 
part  of  this  phase  of  the  work.  They  will  be  established  and 
validated  for  integration  on  an  engine  specific  basis. 

The  next  detection  method  discussed  looks  at  using  a 
commonly  installed  speed  sensor  for  detecting  faults  by 
means  of  measuring  changes  in  the  crankshaft  speed. 

5.  Time  Domain  Speed  Encoder  Analysis 

When  an  engine  misfires,  there  will  be  a  sudden  decrease  in 
engine  speed  immediately  after  the  misfire  event. 
Therefore,  when  a  misfire  event  occurs,  a  torsional  vibration 
is  excited  in  the  crankshaft  and  will  show  up  in  the 
instantaneous  speed  of  the  engine.  The  task  of  this  analysis 
was  to  detect  injector  faults  based  on  changes  in  crankshaft 
speed.  For  the  purpose  of  this  testing,  the  crankshaft  was 
equipped  with  two  speed  encoders,  a  one  pulse  per 
revolution  (1  PPR)  and  a  360  PPR  encoder,  respectively. 
The  one  PPR  encoder  provided  an  index  signal  for  detecting 
each  revolution  of  the  crankshaft.  The  360  PPR  encoder 
enabled  calculation  of  crankshaft  angular  velocity  at  a 
resolution  of  1°  of  revolution. 

In  the  hypothetical  ideal  scenario,  crankshaft  speed  would 
be  constant  with  no  fluctuations  for  a  given  steady  state 
operating  condition.  However  in  practice,  each  firing 
event/impulse  creates  a  sinusoidal  oscillation  of  crankshaft 
speed.  The  blue  line  in  Figure  10  illustrates  this  speed 
oscillation.  The  green,  magenta  and  red  lines  indicate  the 
injector  signals  for  the  respective  cylinders.  Keeping  in 
mind  the  x-axis  represents  720°  degrees  (2  complete 
revolutions  of  the  crankshaft),  the  alternating  yellow  and 
white  regions  corresponds  to  the  firing  order  of  the 
respective  cylinders  #l-5-3-6-2-4. 

Two  additional  comments  are  pertinent  to  interpreting  the 
following  plots: 

1.  The  speed  plot  is  a  composite  mean  of  50  two-revolution 
cycles  of  the  crankshaft.  Averaging  over  multiple  cycles 
over  a  given  steady  state  operating  condition  results  in  the 
averaging  out  of  extraneous  signal  noise.  The  capability  to 
synchronize  multiple  cycles  and  thereby  enable  multi-cycle 
averaging  is  accomplished  by  combining  the  one  PPR 
encoder  channel  as  an  indexing  signal  with  the  injector  one- 
two  channel  (green)  to  identify  a  consistent  crank  angle 
reference  point.  This  reference  point  corresponds  to  the  top 
dead  center  (TDC)  position  of  cylinder  one  at  the  beginning 
of  the  cylinder  one  power  stroke. 


Figure  10.  Sinusoidal  Oscillation  of  Crankshaft  Speed  with 
All  Six  Injectors  Operating 

2.  The  left  vertical  axis  corresponds  to  the  speed  plot  (blue 
line),  and  the  right  vertical  axis  corresponds  to  the  fuel 
injector  signal  voltage.  Signal  levels  greater  than  zero 
correspond  to  the  fuel  injector  injecting  fuel  into  the 
cylinder. 

The  reason  the  injector  signals  are  represented  with  green, 
magenta  and  red  lines  is  due  to  the  instrumentation 
configuration.  There  are  three  pairs  of  injector  voltage 
sensors.  The  injectors  are  grouped  sequentially  (Green: 
injector  one  and  two),  (Magenta:  injector  three  and  four), 
(Red:  injector  five  and  six).  Keep  in  mind,  the  following 
analysis  is  using  the  injector  signals  for  validation  purposes. 
The  motivation  for  this  analysis  assumes  the  fuel  injector 
signals  would  not  be  accessible  as  inputs  for  fielded-real 
world  fault  detection  algorithm. 

The  data  in  Figure  11  shows  the  instantaneous  speed  plots 
over  two  revolutions  of  the  crankshaft  for  a  baseline  (no 
injector  fault)  condition  and  an  injector  number  three  fault 
condition.  At  low  speed  and  load  conditions  (30  lbs-ft  @ 
700  RPM),  the  speed  oscillates  in  a  relatively  sinusoidal 
manner  for  both  the  no-fault  and  fault  conditions.  The 
abnormal  characteristic  of  an  increased  range  in  RPM 
variation  in  the  crankshaft  speed  for  the  faulted  case  (red 
line)  observed  primarily  during  cylinder  two’s  firing  order 
indicates  the  presence  of  an  injector  fault.  Based  on 
observations  across  multiple  low  speed  and  load  setpoint 
conditions,  the  characteristic  ‘spike’  in  speed  oscillation  is 
typically  observed  one  to  two  cylinders  later  in  the  firing 
order. 


6 


Annual  Conference  of  Prognostics  and  Health  Management  Society  2012 


INSTANTANEOUS  SPEED:  LOW  SPEED  (700  RPM);  LOW  LOAD  (30  LBS-FT) 
Mean  of  50  cycle  segments  synchronized  with  respect  to  Cylinder  1  TDC  rising  edge  pulse 


Figure  11.  Low  Speed-Low  Load  Speed  Profile  over  Two 
Revolutions  of  the  Crankshaft 

However  it  was  observed  during  our  investigation  that  the 
potential  to  utilize  this  dynamic  characteristic  for  fault 
detection  purposes  over  the  full  range  of  engine  operation  is 
limited.  This  is  due  to  the  fact  that  the  ‘spike’  is  not  as 
clearly  or  consistently  evident  at  higher  speed  and  load 
conditions.  This  is  discussed  and  shown  in  the  following 
figures. 

The  data  in  Figure  12  and  Figure  13  illustrate  two  examples 
of  the  speed  variation  observed  at  higher  speed  and  load 
conditions  (400  lbs-ft  @  1450  RPM  and  560  lbs-ft  @  2200 
RPM  respectively).  In  Figure  12  the  ‘spike’  is  not  evident.  It 
appears  that  the  relative  drop  in  speed  during  the  crank 
angle  region  corresponding  to  cylinder  six  firing  order  may 
be  indicative  of  the  injector  three  misfire. 


Crank  angle  (Degrees  with  respect  to  start  of  Cylinder  1  window) 


Figure  12.  Speed  Profile  Over  Two  Revolutions  of  the 
Crankshaft  (400  lbs-ft  @  1450  RPM) 


In  Figure  13,  there  is  a  ‘spike’  occurring  during  cylinder 
two’s  firing  order.  However  there  is  also  a  ‘spike’  occurring 
during  cylinder  one’s  firing  order.  These  inconsistent  speed 
variations  require  additional  pre-processing  prior  to 
implementing  automated  fault  detection  algorithms. 


Figure  13.  Speed  Profile  Over  Two  Revolutions  of  the 
Crankshaft  (560  lbs-ft@2200  RPM) 

In  general,  the  correlation  between  fuel  injector  misfire  and 
speed  fluctuation  is  limited  when  using  a  visual  inspection 
of  the  time  domain  speed  plots  based  on  instantaneous 
speed  calculation  of  the  shaft  encoder.  Distinguishing 
injector  misfires  at  higher  speeds  (>700  RPM)  from  ‘noisy’ 
engine  speed  data  is  not  clear-cut  and  will  require  more 
sophisticated  processing  techniques. 

The  instantaneous  speed  based  detection  approach  is 
capable  of  detecting  failure  however  it  is  computationally 
more  complicated  to  implement.  Visual  detection  of  a  given 
injector  fault  is  apparent  at  different  speed  and  load 
conditions.  The  automated  detection  techniques  developed 
in  the  course  of  this  investigation  are  effective  at  detecting 
injection  related  failures  at  low  speeds  and  loads.  At  higher 
speeds,  using  instantaneous  speed  to  distinguish  specific 
injector  faults  is  more  difficult  due  to  the  speed  fluctuations 
resulting  from  torsional  vibration  dynamics  of  the  engine 
crankshaft.  For  these  operating  conditions,  the  pattern 
recognition  algorithms  would  require  more  processing 
power  than  the  other  time-domain  solutions  presented  in  this 
paper. 

For  this  investigation,  this  technique  was  not  considered  for 
further  development  due  to  the  number  of  other  practical 
methods  which  were  found  to  be  suitable  for  embedding 
into  controllers.  However,  instantaneous  speed 
measurements  were  found  to  be  valuable  using  an  order 
domain  approach.  This  technique  enables  detection  of  faults 
using  the  crank  encoder  at  high  speeds  and  loads  as  well  as 
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the  ability  to  identify  the  specific  cylinder  with  the  injector 
fault. 

The  following  sections  will  discuss  the  timer  board  data 
collection  approach  and  the  order  domain  fault  detection 
processing  technique  in  further  detail. 

6.  Order  Domain  Speed  Encoder  Analysis 

This  section  will  discuss  the  method  for  collecting  and 
processing  order  domain  data  for  injector  fault  analysis. 
This  investigation  is  based  upon  work  by  Geveci,  Osburn, 
and  Franchek  (2005). 

6.1.  Data  Collection  Background 

Throughout  the  test  schedule,  the  engine’s  instantaneous 
crankshaft  velocity  was  sampled  using  a  Time  Interval 
Measurement  System  (TIMS).  While  common  A/D 
converters  sample  the  amplitude  of  a  sensor’s  output  at 
discrete  time  intervals,  a  TIMS  records  the  number  of  timer 
board  counts  between  the  output  signal’s  threshold 
crossings.  TIMS  measurement  systems  are  used  in  various 
areas  in  machinery  health  prediction  systems  as  shown  in 
Lebold,  Maynard,  Reichard  (2005),  Lebold  Maynard  (2003) 
and  Maynard,  Lebold,  Groover  ,Trethewey  (2000). 

An  AVL  365C  encoder  system  was  employed  to  produce  an 
output  signal  representative  of  the  crankshaft’s 
instantaneous  velocity.  The  encoder  system’s  periodic 
output  waveform  is  amplified  to  rails  using  threshold  values 
to  produce  a  TTL  logic  signal  having  hi/low  amplitudes 
corresponding  to  the  encoder  wheel’s  optical  pattern  as 
illustrated  in  Figure  14. 


Figure  14.  Conversion  of  sinusoidal  sensor  output  to  square 
wave  pulses. 


The  TTL  logic  levels  are  then  converted  into  a  binary 
representation  to  accurately  define  the  pattern-passing  rate. 
The  exact  transition  time  from  low-to-high  and  high-to-low 
of  an  optical  pattern  passing  is  important  for  processing. 
Therefore  the  fastest  available  analog-to-digital  (A/D) 
converter  would  be  preferable.  However,  these  devices  are 
very  costly  and  are  not  necessary  since  the  replication  of  the 
signal’s  temporal  characteristics  is  not  necessary.  Only 
accurate  capture  of  the  transition  times  are  necessary.  A 
timer/counter  board  is  an  ideal  candidate  for  this  system 
because  it  provides  a  means  of  computing  the  transition  time 
between  an  incoming  logic  signal  with  a  high  degree  of 
resolution. 

The  National  Instruments  PCI-6602  Timer/Counter  Board  is 
a  cost-effective  solution  and  was  used  for  data  collection. 
With  an  80MHz  clock  reference,  pulse-width  measurements 
from  the  PCI-6602  can  be  performed  at  a  resolution  of 
0.0000000125  seconds.  The  data  collection  system  is 
designed  around  the  semi-buffered  period  measurement 
approach  using  direct  memory  access  (DMA)  for  speed.  In 
pulse-width  measurement  applications,  the  board  counts  the 
number  of  ticks  (clock  pulses  of  the  time-base)  that  occur  in 
every  semi-period  of  the  input  waveform.  It  then  stores  this 
value  in  a  data  buffer.  This  results  in  an  array  of  clock 
counts  per  semi-period  of  the  encoder  hub’s  pattern  as 
illustrated  in  Figure  15.  Knowing  the  reference  clock 
frequency  used  during  the  measurement,  the  pulse-width, 
period  or  frequency  of  the  input  waveform  can  be 
computed.  For  this  application,  the  series  of  pulse  widths 
will  be  used  to  determine  the  engine’s  order  content. 


COUNT  1  |[  COOTS  |  COUNT  3  |  COUNT  4  |  COUNT  5  | 


Figure  15.  Zero  crossing  timer  reset  points. 
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6.2.  Order  Domain  Data  Analysis 

Datasets  from  the  TIMS  provided  a  semi-period  resolution 
of  720  PPR.  Before  processing,  all  datasets  were 
normalized  to  begin  at  the  TDC  corresponding  to  the 
beginning  of  cylinder  one’s  power  stroke,  hereafter  referred 
to  as  ‘true  TDC’.  In  this  way,  the  relative  phase  differences 
between  particular  cylinder  firing  events  were  consistent  for 
all  snapshots. 

Encoder  semi-periods  in  the  datasets  were  combined  to 
allow  for  a  period-based  analysis.  Because  encoder  wheels 
typically  have  highs/lows  of  different  widths,  a  period-based 
analysis  eliminates  the  variance  between  semi-periods. 
Shown  below,  Figure  16  illustrates  the  time  history  of  a 
semi-period  signal  (blue)  and  the  constructed  period  based 
signal  (black)  that  was  used  in  the  subsequent  analysis.  It  is 
important  to  note  that  the  mean  of  both  signals  was  removed 
to  illustrate  a  reduction  in  high-frequency  signal  content. 


Example  normalized  encoder  signal  time  history 


Figure  16.  The  Semi-Period  Signal  (Blue)  Recorded  in 
Snapshot  492  During  Test  #25  and  the  Constructed  Period 
Based  Signal  (Black) 

The  period  based  signals  were  then  demodulated  to 
calculate  the  amplitude  of  the  crankshaft’s  torsional 
vibration  in  degrees  in  the  rotational  domain.  To  reduce  the 
effects  of  events  uncorrelated  to  the  shaft’s  rotation  in  the 
demodulated  signal,  synchronous  averaging  was  also 
applied  in  the  rotational  domain.  The  amplitude  and  phases 
of  the  demodulated,  averaged  signal  were  then  analyzed  in 
the  order  domain. 

Figure  17  through  Figure  19  illustrate  order  spectra 
comparing  data  collected  during  baseline  and  injector  fault 
tests  at  the  discrete  low,  medium  and  high  engine  speeds, 
respectively. 


Order  spectrum:  diesel  engine  data 


Shaft  order 

Figure  17.  Baseline  (Snapshot  492)  and  Cylinder  Three 
Fault  (Snapshot  487)  Order  Spectra  for  Fow  Speed,  Fow 
Toad  Engine  Operating  Conditions 

An  order  domain  analysis  of  the  data  collected  during  the 
injector  fault  test  schedule  suggests  that  cylinder  faults  can 
be  identified  by  monitoring  the  amplitude  of  the 
crankshaft’s  orders  lower  than  two.  When  an  injector  fault 
was  introduced  into  the  test  system  at  low,  medium,  and 
high  load/speed  engine  conditions,  noticeable  modulation  of 
the  crankshaft  occurred  between  orders  one-half  and  two. 


Order  spectrum:  diesel  engine  data 
1450  RPM,  33.7%,  36  ft-lbs 


Shaft  order 


Figure  18.  Baseline  (Snapshot  446)  and  Cylinder  Three 
Fault  (Snapshot  443)  Order  Spectra  for  Medium  Speed, 
Medium  Toad  Engine  Operating  Conditions 

Using  these  same  datasets,  an  injector  fault  threshold 
amplitude  for  the  orders  between  one-half  and  two  was 
calculated.  When  the  amplitude  of  any  of  the  orders 
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between  one-half  and  two  rises  above  this  threshold  value 
and  enters  the  “fault  detection  region”,  a  cylinder  fault  is 
diagnosed.  The  injector  fault  detection  region  above  the 
empirically  determined  0.003  degrees  threshold  is  illustrated 
in  Figure  17  through  Figure  19. 


Order  spectrum:  diesel  engine  data 
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Shaft  order 

Figure  19.  Baseline  (Snapshot  485)  and  Cylinder  Three 
Fault  (Snapshot  479)  Order  Spectra  for  High  Speed,  High 
Load  Engine  Operating  Conditions 

The  effect  of  an  injector  fault  on  the  amplitude  of  the  half 
order  can  be  seen  in  Figure  17  through  Figure  19. 

Figure  19  testing  results  suggest  that  the  magnitude  of  the 
torsional  modulation  caused  by  an  injector  fault  is 
dependent  on  the  engine’s  operating  conditions.  At  low 
speed,  low  load  conditions  the  injector  fault  causes  little 
torsional  modulation  to  the  shaft’s  steady  state  rotation  at 
the  half  order.  At  high  speed,  high  load  conditions, 
however,  the  torsional  modulation  at  the  one-half  order  is 
very  pronounced. 

After  an  injector  fault  is  diagnosed,  it  is  necessary  to 
identify  the  cylinder  that  contains  the  injector  fault.  This, 
too,  was  accomplished  in  the  order  domain.  A  shaft  order 
can  be  defined  to  describe  the  frequency  of  a  particular 
synchronous  event  that  occurs  at  each  shaft  revolution. 
From  this  definition,  one  can  conclude  that  a  shaft’s  one 
half  order  contains  information  about  events  that  occur  once 
every  two  shaft  revolutions.  For  a  six  cylinder,  four- stroke 
engine  this  order  specifically  characterizes  the  frequency  of 
a  crankshaft’s  torsional  loads  that  result  from  an  injector 
fault.  These  loads  are  caused  by  the  energy  loss  due  to  lack 
of  cylinder  combustion,  which  in  turn  results  in  a 
momentary  decrease  in  crankshaft  velocity  relative  to  the 
no-fault  steady  state  crankshaft  velocity  profile. 


To  identify  an  injector  fault,  the  phase  of  the  half  order  was 
analyzed  for  each  dataset.  Because  all  datasets  were 
normalized  to  begin  with  the  first  true  TDC  of  the  first 
cylinder  in  the  engine’s  firing  sequence,  the  relative  phase 
of  cylinder  firing  events  was  maintained  in  this  analysis. 
After  the  order  spectrum  was  calculated  for  each  dataset,  the 
phase  of  the  one-half  order  was  extracted.  Figure  20 
illustrates  the  phase  of  each  dataset’s  half  order  as  injector 
signals  to  various  cylinders  were  deliberately  attenuated. 
The  phase  clustering  of  the  one-half  order  provides  a  means 
to  identify  a  particular  faulty  injector  signal. 

An  analysis  of  Figure  20  reveals  that  the  one-half  order 
phases  are  clustered  in  groups  that  are  separated  by  roughly 
120°.  The  physical  significance  of  this  separation  angle  can 
be  traced  to  the  crank  rotation  angle  between  cylinder  firing 
events  that  is  120°  for  a  four-stroke,  six  cylinder  engine. 


Phase  clustering  for  diesel  engine  cylinder  diagnostics 
Crankshaft's  1/2  order  analyzed 


Figure  20.  Phase  Clustering  of  One-Half  Order 

Order  domain  inspection  of  data  collected  by  the 
crankshaft’s  existing  encoder  system  provides  valuable 
information  that  can  be  used  for  injector  fault  diagnostics. 
By  monitoring  the  amplitude  of  the  crankshaft’s  low  order 
content  relative  to  the  suggested  fault  detection  window, 
injector  health  can  be  classified.  If  a  faulty  injector  signal  is 
diagnosed,  the  phase  of  the  one-half  order  can  be  compared 
to  the  clustering  in  Figure  20  to  determine  the  cylinder  that 
is  most  likely  faulted. 

7.  Conclusion 

This  paper  described  four  approaches  for  detecting  fuel 
injector  misfire  events;  an  injector  signal  based  approach,  a 
vibration  signal  approach,  a  crankshaft  speed  based 
approach,  and  an  order  domain  approach.  Relative  to  the 
crankshaft  speed  based  approach,  the  injector  signal  and 
vibration  time  domain  approach  both  showed  the  greatest 
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potential  for  on-board  vehicle  implementation.  These 
approaches  demonstrated  the  capability  of  detecting 
individual  injector  misfires  on  the  seven  liter  diesel  engine 
over  the  (30  lbs-ft  @  700  RPM)  to  (560  lbs-ft  @2100 
RPM)  operating  range  when  tested  on  an  engine 
dynamometer.  Furthermore,  each  approach  equally 
demonstrated  the  ability  to  correctly  identify  the  specific 
cylinder  location  of  the  faulted  injector. 

Both  approaches  require  a  TDC/one  PPR  or  360PPR 
encoder  signal  for  synchronization  or  the  camshaft  speed 
encoder  signal,  along  with  the  basic  hardware  and  software 
for  algorithm  implementation. 

The  injector  signal  based  approach  does  require  access  to 
the  signal  taps  for  the  six  respective  injectors  while  the 
vibration  based  approach  requires  one  accelerometer 
mounted  to  the  top  of  the  engine.  Due  to  the  nature  of 
vibration  monitoring,  this  approach  should  also  work  for 
detecting  any  combustion  related  faults  but  must  be  first 
validated  on  an  actual  platform  to  ensure  it  is  immune  to 
externally  induced  vibrations.  Among  other  factors,  the 
extent  to  which  fuel  line  clogs  or  worn  injector  actuators  are 
relevant  faults  may  dictate  the  reliability  of  this  approach  in 
a  fielded  environment. 

The  order  domain  approach  showed  great  potential  for  on¬ 
board  vehicle  implementation.  The  approach  demonstrated 
the  capability  of  detecting  individual  injector  misfires  on  the 
seven  liter  diesel  engine  over  the  (30  lbs-ft  @  700  RPM)  to 
(560  lbs-ft  @2100  RPM)  operating  range  when  tested  on  an 
engine  dynamometer.  Additionally,  this  approach 
demonstrated  the  ability  to  correctly  identify  the  specific 
cylinder  location  of  the  faulted  injector. 
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